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1.  Introduction 


Analysis  and  predictions  of  projectile  trajectories  has  a  long  history.  Notable  authors  are 
Galileo,  Bernoulli,  and  Euler  who  found  mathematical  solutions  to  this  problem  which  gave 
important  advances  to  trajectory  predictions.  Newton  has  been  credited  for  the  quadratic  law  of 
resistance  characterizing  the  aerodynamic  drag  of  a  body.  Further  study  has  shown  that  a  more 
sophisticated  characterization  of  drag  is  required  to  predict  accurate  trajectories.  Weinacht 
et  al.  (1)  give  a  nice  historical  perspective  of  previous  work  followed  by  a  detailed  analysis  of 
trajectories  based  on  power-law  drag  curves.  Increased  sophistication  has  also  lead  to  some 
unfortunate  consequences.  Details  are  often  required  that  may  not  always  be  relevant  to  the 
problem  in  question  and  thus  over-burdening  the  user  or  designer.  For  such  cases,  simplified 
analysis  can  provide  accurate  results  with  a  minimum  of  relevant  input. 

High-velocity,  direct-fire  munitions  can  be  studied  using  point  mass  trajectory  equations  coupled 
a  with  power-law  drag  relation  (1-4).  The  study  (i),  suggested  by  the  experimental  study  of 
Celmins  (5),  presented  analysis  and  predictions  of  flat- fire  projectiles  with  small  yaw  and  high 
velocities  using  a  drag  law  of  the  form: 


C  =C 

*-0  v„ 


vVy 


(1) 


The  governing  equations  assumed  two-dimensional  (2-D)  motions  with  parameterized  drag 
curves  depicting  drag  for  projectile  velocities  that  result  in  high  Mach  numbers.  The  analysis 
was  done  by  first  neglecting  gravity  and  finding  closed-form  solutions  to  the  zero  gravity 
velocity  profiles.  Then  motion  under  the  influence  of  gravity  was  modeled  using  drag  forces 
(equation  1)  subjected  to  these  zero  gravity  velocities. 

A  continuation  of  this  work  is  given  here  but  this  analysis  retains  gravity  such  that  drag  is  always 
influenced  by  gravity  while  using  the  same  power-law  drag  formulation  found  in  Weinacht  et  al. 
(1).  The  vertical  and  horizontal  momentum  equations  are  shown  to  be  equivalent  to  a  single 
third-order  differential  equation  (DE)  in  which  the  independent  variable  is  the  down-range 
distance.  Making  the  assumption  that  the  explicit  slope  of  the  trajectory  DE  can  be  written  as  a 
free  parameter  allows  the  governing  third-order  DE  to  be  integrated  giving  simple  closed-form 
solutions.  When  the  launch  angle  is  small,  these  solutions  become  the  solutions  found  in 
Weinacht  et  al.  (1).  A  value  for  the  free  parameter  is  given  and  is  shown  to  produce  very 
accurate  closed-form  solutions  for  launch  angles  up  to  30°.  These  solutions  yield  easy-to-use 
results  that  broaden  those  given  in  Weinacht  et  al.  (1)  for  fiat-fire  cases.  The  theory  and  model 
discussions  are  presented  in  the  following  sections. 
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2.  Equations  of  Motion 


Assuming  the  2-D  theory  found  in  Weinaeht  et  al.  (1)  gives  the  governing  equation; 

dt 


(Y  V 

x  =  -},V  ^  i 


fY  Y 

y  =  -},V  ^  y-g 


V'  =  x"  +  f 


Using  X  as  the  independent  variable  leads  to  y  =  y'  x  and  y  =  y"  x^  +  y'  x  ,  so  equation  2 
beeomes: 


y  ^ 

X 


m  2gX 

y  =^^  =  - 

X 


2AgV  ^ 


Noting  that  V  =  x  Jl  +  y'"  =  — gives  the  following  expressions  and  initial  conditions: 


_  _  n+z 

.  zvi-vlJ-ZU  ,  „  .1  .  a-,  .1 

y  = - s — ^ — su —  yl  =  0  ’  y  n  =  ’  y  L  =  ■ 


[y"+i]' 


Vocos"4) 


x=  — ^  ,  y  =  xy'  ,x|  =0  andV'l  . 

y  y 


For  eonvenienee,  two  length  seales  are  defined  asy  =  Y g/V'|^^  and  x  =  X  Vq/V'I^  making  the  last 
equation  in  equation  3  take  the  form 


2 
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(5) 


Y"'  =  —  - 


.  - n+2  I 

2(-Y-rT  _  I  V.V'Ltan^,  __ 

n-1  ’  ^  n  ^  n  _  ’  ^  n 


gY' 

VoV'L 


2  A 

+1 


Vq  cos  (j) 


Please  note  that  equation  5  is  equivalent  to  equation  1  even  when  =  >b(y) ,  but  the  presentation 
given  here  takes  X  =  eonstant . 

The  analysis  given  in  Weinaeht  et  al.  (1)  foeused  on  direet  fire,  thus  assuming  the  elevation  angle 

n2 

gY' 


(j)  is  small.  This  suggests  a  very  good  approximation  is  to  let 


V.Vl 


« 1 ,  whieh  leads  to 


results  nearly  identieal  to  those  given  in  Weinaeht  et  al.  (1)  if  (j)  is  small.  However,  in  the 
interest  of  using  elevation  angles  (j)  that  are  not  neeessarily  small,  it  further  suggests  treating 


Q  = 


gY' 

VoV'L_ 


2  A 

+1 


n-1 


as  a  free  parameter  whieh  allows  the  following  simple  elosed-form 


solution  to  equation  5  satisfying  the  initial  eonditions. 


Y  =  -Q 


r 


Qeos""  > 


Xeos”  (j) 
(n-2) 


Xn  ^ 
y  Q  eos”  (j)  y 


2 11-2  ^ 


-1 


+  - 


V 


2(n  -l)(n  -2) 


+  - 


tanc^VoV'lgX 


(6) 


Note  the  last  expression  beeomes  identieal  to  the  results  presented  in  Weinaeht  et  al.  (1)  when 
(j)  «  1 ,  thus  making  eos  (j)  — >  1  and  Q  ^  1 .  Sinee  the  remainder  of  this  report  stems  from 

manipulations  of  equation  6,  it  eauses  all  following  results  given  here  to  have  similar  identieal 
behavior  to  Weinaeht  et  al.  (1)  whenever  the  launeh  angle  beeomes  small.  This  report  will  show 
equation  6  gives  good  solutions  to  the  parameterized  equation  5  for  elevation  angles  0  <  (j)  <  30° 

by  using  Q  =  eos' "  ^  .  This  ehoiee  for  Q  was  obtained  using  trial  and  error  while  eomparing 

numerieal  solutions  of  equation  5  to  those  of  equation  6,  starting  with  Q  =  1 .  The  authors  have 
judged  the  given  value  for  Q  is  adequate  for  0  <  (j)  <  30°.  The  following  presentation  primarily 
foeuses  on  (j)  =  30°  sinee  these  eases  were  found  to  display  the  weakest  agreement  of  all  launeh 
angles  in  the  range  0  <  (j)  <  30°.  However,  a  few  exeursions  with  (j)  =  45°  are  diseussed  to 

indieate  the  limitations  of  Q  =  eos' "  —  for  (j)  >  30°. 
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For  completeness,  the  limiting  cases  for  equation  6  are  given  in  the  following  expressions: 


Y- 


tan(^VoV'|„X 


f  2xA 


Q  +  2X-Qe' 
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Q"  log 


4  eos  (j) 
X"  2X 
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X 
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log 
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2X 
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(V) 


The  ehoiees  of  projeetiles  are  the  same  as  those  used  in  Weinaeht  et  al.  (1)  and  the  pertinent  data 
are  taken  from  table  1 . 


Table  1.  Projectile  characteristics. 


Projectile  Type 

Muzzle  Velocity 

Vo 

(m/s) 

Muzzle  Retardatiou 

([m/s]/km 

V'o 

Vo 

(1/km) 

M829A1 

1580 

68 

0.043 

M865PIP 

1700 

343 

0.202 

M830 

1140 

273 

0.239 

M830A1 

1410 

209 

0.148 

Comparative  examples  illustrating  the  validity  of  the  parameterized  solution  equation  6  to  the 
numerieal  solution  of  equation  5  are  given  in  figure  1 .  The  agreement  between  the  two  types  of 
solutions  is  very  strong  except  where  the  trajectory  slope  is  steep  for  the  M865PIP  projeetile. 

The  slope  of  the  trajectories  shown  in  figure  1  is  given  in  figure  2,  where  it  is  evident  that  the 
analytic  solutions  agree  with  the  numerieal  solutions.  However,  the  partieular  case,  M865PIP, 
has  weaker  agreement  but  this  can  be  attributed  to  the  mueh  steeper  slope  and  the  fact  that 
derivatives  of  an  approximation  usually  don’t  agree  as  well  as  the  given  approximation. 
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Figure  1.  Comparison  of  analytic  and  numeric  trajectories. 


Figure  2.  Comparison  of  analytic  and  numeric  trajectory  slopes. 

Velocity  components  are  readily  obtained  from  equation  6  and  limiting  values  from  equation  7 
produce  the  following  expressions; 
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gQlog 


f  2XA 

1-e^ 

V  J 
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1- 


2X 
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2V  Vl 
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Figure  3  shows  Mach  number  M  profiles  calculated  from  the  last  set  of  equations  where 
comparisons  are  made  to  numerieal  solutions  generated  from  equation  5. 


Figure  3.  Comparison  of  numerical  and  analytic  Mach  number  profiles. 


(8) 
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The  last  set  of  plots  also  show  that  the  weakest  agreement  of  Y  and 


dY 


figures  1  and  2, 


occurs  when  the  Mach  number  M<1 .  Since  the  power- law  drag  formulation,  equation  2,  is 
usually  considered  for  M>1,  it  suggests  that  the  given  results  for  (j)  <  30°  are  more  applicable  for 
ranges  where  M>1 .  Examples  of  the  types  of  errors  that  occur  when  (j)  >  30°  are  shown  in 
figure  4,  where  (j)  =  45°. 


Trajectory  <!>  -  45 
Red  n=.75  Blue  n=.25  Black  n=0 


50 
45 
40 
35 
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Y  25 
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0.0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6  1.8 


X 


Figure  4.  Comparison  of  analytic  and  numeric  trajectories  for  cp  =  45. 

These  results  may  be  interpreted  as  an  indication  that  0  <  (j)  <  30°  is  an  over-restriction  for  the 
ranges  where  M>1 . 


3.  Gravity  Drop 


Analytical  results  of  the  previous  work  (7)  regarding  gravity  drop  were  obtained  by  considering 
dimensional  y  terms  that  are  dependent  on  g.  Applying  this  last  statement  to  the  analysis  given 
here  simply  means  to  set  the  term  proportional  to  g~^  in  equations  6  and  7  to  zero  to  get  scaled 
versions  of  gravity  drop  “Yg_drop.”  Furthermore,  expanding  this  expression  in  a  Taylor  series  and 
converting  the  series  to  a  rational  function,  thus  increasing  the  range  of  validity,  yields  the 
following  expression; 
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(9) 


cos”  >  Q  X"  (2  (3  n  +  l)x  - 1 5  cos  (^) 

~  3  ( (n  +  2)X  (X  -  4  Q  cos”  ^  )+ 1 0  Q"  cos'"  c^)  ’ 

Figure  5  has  examples  of  gravity  drop  for  the  M865PIP  projeetile  at  various  values  of  the  drag- 
eurve  parameter  n.  Evidently,  is  well  approximated  by  equation  9. 


M865PIP  <t)  =  30° 


Figure  5.  Comparison  of  numerical  and  analytic  gravity  drop  for  the  M865PIP. 


4.  Gun  Elevation  and  Target  Position 


Predicting  the  gun  elevation  angle  (j)  to  aequire  a  target  at  Y  =  0 ,  for  a  given  range  position  X , 
is  a  useful  flight  parameter  (7)  that  can  be  found  from  equation  6.  Using  the  result  for  Yg_^gp 

given  in  equation  9  shows  the  elevation  angle  satisfies 

VoV'lptanct)  ^  -eos”^>QX(2(3n  +  l)X-15eos^) 

g  ~3((n  +  2)x(X-4Q  cos”  (t))+10Q^  cos^ 

Solving  this  expression  for  (j)  as  large  as  30°  generally  will  require  numerieal  techniques.  But 
limiting  analysis  to  small  angles,  (j) «  1  ^  eos  (j)  — >  1  and  Q  — >  1 ,  allows  equation  10  to  give  direet 
solutions  for  small  (j).  These  values  are  the  same  as  those  found  in  Weinaeht  et  al.  (7)  and  are 
repeated  here  for  eompleteness  in  figure  6. 
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Figure  6.  Scaled  gun  elevations  as  a  function  of  scaled  range. 


5.  Trajectory  Rate  of  Change  Due  to  Muzzle  Velocity 


The  rate  of  vertical  change  of  trajectory  due  to  muzzle  velocity  fj,  can  change  significantly  with 
range.  An  analytic  expression  for  this  comes  directly  from  equation  6  and  the  dimensional  rate 

of  change  scaled  as  — — =  Vg  is  shown  in  the  following  expressions; 
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d  Y  _  cos^"~^ 
dVg  n(n-2) 
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1-n- 
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2n-2 


n-2 
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1-n- 
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Qcos”  (j) 


+  n- 
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Qcos”  (j) 


QX 


2X  A 

1-e^ 

V  y 


2cos^  (j) 


n  ^  0 


QX 


log 


1-- 


2X 

Qcos^  (j) 


n  ^  2. 


(11) 


To  validate  equation  11,  the  derivative  of  equation  5  with  respect  to  Tg  is  numerically  integrated. 

Results  given  in  figure  7  compare  the  parameterized  analytic  solution  to  the  nonparametric 
numerical  solution. 
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Figure  7.  Comparison  of  scaled  numerical  and  analytic  dy/dVo- 

Figure  7  shows  equation  1 1  is  a  good  approximation  to  the  rate  of  ehange  of  the  trajeetory 
caused  by  varying  the  muzzle  velocity  Vg . 


6.  Trajectory  Rate  of  Change  Due  to  Muzzle  Retardation 


The  rate  of  vertical  change  in  the  trajectory  as  a  function  of  the  drag  coefficient  or  muzzle 
retardation  V'l^  can  also  be  an  important  influence  on  trajectory  that  depends  on  range.  Again, 

equation  6  is  used  to  find  a  scaled  version  of  this  rate  of  change. 


dY 

dV'L 


Q2  „^„2n-2  1 

COS  q) 

n  -2 


f 

1 

V 


nX 

Qcos" 

n-1 


+  - 


X 


Qcos”  (j) 


n-2 

nX 

Qcos"  (t)y 


+  1 


-Q 


2  cos  (j) 


2X 


(X-Q)e‘^  +X  +  Q 


-2QMog 


1-- 


X 


QX 


n  ^  0 
X 

Qcosct) 


-2 


QX- 


QcOSCt)  y 
cos^  (j) 


r 


X 


Qcos  (j) 


coscj) 


--1 


X 
Qcosct) 


,n  ^  1 


-1 


log 


1-- 


2X 


Qcos  (j) 


n  ^  2. 


(12) 


10 


Comparing  these  expressions  to  the  numerieal  integration  of  the  derivative  of  equation  5  with 
respect  to  V'l^  is  shown  in  figure  8.  The  case  n  =  0.25  shows  the  analytic  solution  degrades  with 

increasing  range.  Once  again,  this  is  not  too  surprising  since  the  derivative  of  a  well-converged 
approximate  function  suffers  with  less  convergence. 


Figure  8.  Comparison  of  scaled  numerical  and  analytic  d  y/d  . 


7.  Time  as  the  Independent  Variable 


The  presentation  given  here  has  considered  x  as  the  independent  variable  but  all  of  the  previous 
results  could  be  given  in  terms  of  the  time  t  if  desired.  Integrating  the  first  equation  in  equation  2 
leads  to  an  expression  of  t  as  a  function  of  the  scaled  range  X  given  by 
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Validation  of  equation  13  is  displayed  in  figure  9,  where  it  is  evident  that  this  expression 
supplies  an  accurate  method  to  covert  the  scaled  spatial  independent  variable  X  to  the  scaled 
temporal  variable  t  .  Some  error  is  introduced  for  the  M865PIP  projectile  which  can  be 
attributed  to  the  rapidly  decreasing  slope  for  X  >  1  (see  figure  2).  This  may  not  be  a  serious 
defect  due  to  the  relative  small  Mach  number  for  X  in  this  range. 


Figure  9.  Comparison  of  scaled  time  as  a  function  of  scaled  range. 


8.  Low  Velocity  Projectile 


All  of  the  investigations  in  this  report  have  thus  far  focused  on  projectiles  with  high  Mach 
number  launch  velocities.  Therefore,  suspicions  may  arise  when  applying  this  theory  to 
projectiles  with  low  launch  velocities.  To  address  this  issue  implies  considering  a  sample  of  a 
relatively  low  launch  velocity  projectile.  One  such  projectile  is  Scorpion,  which  has  a  typical 
launch  velocity  of  Vo  =  76  m/s  and  a  muzzle  retardation  of  V'j^  =  104([m/s]/km) .  Figures  10  and 
1 1  have  examples  of  subjecting  Scorpion  to  the  theory  of  this  report  using  the  usual  value  of 
n  =  0  for  subsonic  launch  velocities. 

Evidently,  the  parametric  results  for  low  velocities  agree  as  well  as  they  do  for  high  Mach 
number  projectiles  with  launch  angles  0  <  (()  <  30°.  When  (j)  =  45°,  the  agreement  has  degraded 
but  this  agreement  may  still  be  strong  enough  to  suggest  that  (j)  <  30°  is  somewhat  over- 
restrictive. 
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Scorpion  n=0 


X 


Figure  10.  Comparison  of  trajectory  calculations  for  Scorpion. 


Scorpion  n=0 


X 


Figure  11.  Comparison  of  trajectory  slope  calculations  for  Scorpion. 
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9.  Summary 


This  report  extends  the  investigation  found  in  Weinaeht  et  al.  (1)  whieh  eonsidered  direct-fire, 
high-velocity  projectiles  subjected  to  power-law  drag  curves.  The  extension  here  continues  to 
assume  power- law  drag  curves  but  now  the  launch  angle  (j)  takes  on  values  as  large  as  (j)  =  30°. 

In  Weinaeht  et  al.  (1),  gravity  is  treated  as  a  small  perturbation  when  projectile  velocities  are 
high.  However,  in  this  report,  the  equations  of  motion  are  written  as  a  third-order  nonlinear 
differential  equation,  DE,  where  range  is  the  independent  variable.  A  simple  analytic  solution  to 
the  DE  is  found  by  assuming  a  factor  in  this  DE  can  be  treated  as  a  fixed  parameter  Q  .  A  value 

for  Q  =  cos^  “  ^  is  shown  to  produce  valid  solutions  for  launch  angles  0  <  (j)  <  30°  and,  in  the 

limit  of  (j)  ^  0 ,  yields  solutions  identical  to  those  given  in  Weinaeht  et  al.  (1).  Experience  has 
shown  the  validity  of  the  parameterized  solutions  increases  as  (j)  becomes  smaller  and,  in  fact,  if 
Q  =  1  for  0  <  (j)  <  20°,  produces  accurate  solutions.  Examples  where  (j)  =  45°  characterize  the 
errors  and  trajectory  degradation  that  occurs  when  the  launch  angle  is  greater  than  30°,  which 
may  further  indicate  that  0  <  (j)  <  30°  may  be  too  restrictive.  Replacing  the  spatial  dependency  to 
temporal  dependency  is  easily  accomplished  by  application  of  equation  13.  Results  given  here 
can  be  used  in  engineering  applications  that  are  similar  those  discussed  in  Weinaeht  et  al.  (1),  but 
now  the  launch  angle  (j)  can  be  as  large  as  (j)  =  30°. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


D 

g 

m 

M 


n 

S 


t 

V 


X 

y 


p 


Drag  coefficient 
Launeh  drag  coefficient 

Reference  diameter 
Gravitational  acceleration 

Projectile  mass 
Maeh  number 

Exponent  of  drag  power  law 
Reference  area  S  =  7t:D^/4 
Time 

Total  velocity 
Muzzle  velocity 

Dimension  down-range  distanee 
Dimension  vertical  distance 

Muzzle  retardation 
Initial  gun  elevation  angle 
Atmospheric  density 
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